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ABSTRACT 

We present a variability analysis of the early-release first quarter of data publicly released by the 
Kepler project. Using the stellar parameters from the Kepler Input Catalog, we have separated the 
sample into 129,000 dwarfs and 17,000 giants, and further sub-divided the luminosity classes into 
temperature bins corresponding approximately to the spectral classes A, F, G, K, and M. Utilizing 
the inherent sampling and time baseline of the public dataset (30 minute sampling and 33.5 day 
baseline), we have explored the variability of the stellar sample. The overall variability rate of the 
dwarfs is 25% for the entire sample, but can reach 100% for the brightest groups of stars in the 
sample. G-dwarfs are found to be the most stable with a dispersion floor of cr ~ 0.04 mmag. At 
the precision of Kepler, > 95% of the giant stars are variable with a noise floor of ^ 0.1 mmag, 
0.3 mmag, and 10 mmag for the G-giants, K-giants, and M-giants, respectively. The photometric 
dispersion of the giants is consistent with acoustic variations of the photosphere; the photometrically- 
derived predicted radial velocity distribution for the K-giants is in agreement with the measured radial 
velocity distribution. We have also briefly explored the variability fraction as a function of dataset 
baseline (1-33 days), at the native 30-minute sampling of the public Kepler data. To within the 
limitations of the data, we find that the overall variability fractions increase as the dataset baseline is 
increased from 1 day to 33 days, in particular for the most variable stars. The lower mass M-dwarf, 
K-dwarf, G-dwarf stars increase their variability more significantly than the higher mass F-dwarf and 
A-dwarf stars as the time-baseline is increased, indicating that the variability of the lower mass stars 
is mostly characterized by timescales of weeks while the variability of the higher mass stars is mostly 
characterized by timescales of days. A study of the distribution of the variability as a function of 
galactic latitude suggests sources closer to the galactic plane are more variable. This may be the 
result of sampling differing populations (i.e., ages) as a function of latitude or may be the result 
of higher background contamination that is inflating the variability fractions at lower latitudes. A 
comparison of the M dwarf statistics to the variability of 29 known bright M dwarfs indicates that 
the M dwarfs are primarily variable on timescales of weeks or longer presumably dominated by spots 
and binarity. But on shorter timescales of hours which are relevant for planetary transit detection, 
the stars are significantly less variable, with ~ 80% having 12-hour dispersions of 0.5 mmag or less. 
Subject headings: stars: variable, stars: statistics 



1. INTRODUCTION 

Stars have been known for a long time to vary in bright- 
ness, and photometric studies over the past centuries 
have revealed many classes of stars exhibiting a variety 
of variability (Pickering^ ^188ip . With interest in stel- 
lar variability growing tremendously in the last decade 
as ground-based and space-based surveys for exoplanets 
have gained momentum, understanding the stellar pho- 
tometric variability is even more crucial. 

Sources of stellar variability i nclude pulsatio ns, bina- 
rity, rotation, and activity (e.g., |Eyer fc Mowl avi 2008]). 
Having a large sample of uniformly observed stars is vital 
in the categorization and characterization of the variabil- 
ity which can inform us about the stars themselves, their 
companions and companion rates, and their evolution. 
The fractions of stars that are found to be variable is de- 
pendent upon the sample of stars studied, the precision of 
the survey, the range of magnitudes over which the preci- 
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sion is m atched, and the time duration of the survey (e.g.. 



Eyer fc Mowlavin2008i ,Howell. 2008) . For example, Hip- 



parcos (with mmag precision and a completeness limit 
near V=8 mag) found 10% of the stars in the s ample to 
be variable ( |Eyer fc Grenon|1997|[Eyer fc Mowl avi 2008), 
but the variability traction depended upon both the stel- 
lar brightness and the stellar type. Similar population 
and precision dependent results have been found by sur- 
vey programs intended for other purposes such as mi- 



crolensing studies and t ransit surv eys (e.K.,OGLE (Woz 
niak fc Szy manski'1 998i HATN et (|Hartman et al.'2004'). 



and WASPO ( Kane et al.||200 5|)) asjwell as from general 



variabi lity programs (e.g., B SVS (Everett et al. 2002)^^ 
FSVS (IHuber , Everett, fc HoweU 20061, and ASAS LFoA 
jmanskipUMI ): 



As the surveys have become more sensitive, the frac- 
tion of stars observed to vary has been found to increase 
in a form which can be described by a power-law distribu- 
tion directly proportional to the quality of the photomet- 
ric precision powell)2008[ ). This is a result of the current 
"best" survey precisions, time samplings, and survey du- 
rations probing ever deeper into the variability of stars 
but generally not reaching the astrophysical variability 
floor. 
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S paced -based mi ssions such as MOST (Matthews et 
ar||1999|), CoRoT (lA^vergne et aL||2009[ ), and Kepler 
iBorucki et al. 201(jp talte advantage ot the controlled 
environment m space to achieve the best possible pre- 
cision for the telescope ~ increasing the precision of the 
photometry and allowing us to explore the limits of stel- 
lar variability. The Kepler mission, with its large Im 
aperture and huge focal plane (~ 100'-'°), is obtaining 
sub-millimagnitude precision (30 minute integration) and 
micro-magnitude precision (6 hour timescale) for thou- 
sands of stars and has the potential to expand our knowl- 
edge of the limits of stellar variability. 

Kepler was launched in March 2009 and began science 
operations in May 2009. Like CoRoT, Kepler does not 
study all stars within its field-of-view, but rather Ke pler 
monitors a specific set of ^ 150, 000 target stars ( jBatalha| 
et al.|2010 1 . Early work on the variability of stars in the 
Kepler dataset has been performed; these works have 
concentrated on the dwarf stars, periodicity , and flares 



P JWalkowicz et al.|[20To| . 

eKepier project released tc 



( |Basri et al ||2010a|] 

In June 2010, the kepier project released to the public 
the first major time series data product for the majority 
of the targets. We present a discussion of the dataset 
(§ 2.1) and how it it is divided into spectral and lumi- 
nosity classes (§ 2.2). We primarily discuss the stellar 



variability of the sample on the time scale of the dataset 
(33 days) and at the sampling rate of the data (30 min- 
utes); we do explore briefly the variability as a function 
of the time baseline from 1-33 days. Discussions of 



the stellar photometric dispersions (§ 3.1) and the vari- 
ability fractions (§ 3.2) for the dataset as a whole (30- 
minute sampling, 33.5 day baseline) are presented. The 
variability study is extended by exploring the source of 
the variability in the giant stars, the time-dependency of 
the variability fractions, and the vari abili ty fraction as 
a function of galactic distribution (§ |3.3[ ). Finally, we 
explore in more detail the v ariability of the lower mass 
main sequence stars (§ 3.4). Studies and characteriza- 



tion of stellar variability not only provide insight into 
the nature of stars themselves, but also help inform our 
statistical understanding of the detection of transiting 
exoplanets in the presence of stellar "noise" . 

2. KEPLER PUBLIC DATA 

2.1. Quarter 1 and Characterization 

The Kepler project publicly released light curve data 
for all targets observed in the first two "quarters" of ob- 
serving (QO and Ql) and for targets listed by the Kepler 
project as "dropped" from observation in quarters QO, 
Ql, and Q3. We have chosen to utilize only the Ql data 
for this study, as these data represent the most com- 
plete and most uniform set of Kepler data available to 
the public. The Ql data mark the beginning of science 
operations and span approximately 33.5 days from the 
end of QO (13 May 2009) to first spacecraft roll (15 June 
2009^1^ We have also chosen to use only the 30 minute ca- 
dence data (and not the 1 minute cadence data) to main- 
tain the uniformity and continuity of the sample. The 
data are available through the Kepler mission archive at 
MASlQand also through the NASA Star and Exoplanet 

^ Data released to the public June 2010. 
® http:/ /archive. stsci.edu/kepler 



Database (NStED)Q 

In addition to providing access to the light curve data 
themselves, NStED calculates a standard set of statistics 
for each light curve as a whole (33 day baseline at 30 
minute sampling) including a median value, a median of 
the uncertainties, a dispersion about the median value, 
and a reduced chi-square assuming a constant (median) 
value. The statistics are provided as part of the header 
information in the NStED ASCII versions of the pub- 
lic FITS files, and are also searchable and download- 
able as part of the NStED data query service. These 
statistics are calculated on the data corrected by the Ke- 
pler project for "instrumental effects" (ap-Cor r_flux). As 
ment ioned in the Kepler Data Release Notes (van Cleve 
20101, the Kepler project is in the early development 



stages of the data processing pipeline, which is primarily 
intended to find exoplanetary transits. The pipeline may 
not perfectly preserve general stellar variability with am- 
plitudes comparable to or smaller than the instrumental 
systematics on long timescales. 

The Kepler project warns that trends in the data com- 
parable to the length of the time-series data (~ 20 — 30 
days in the case of the Ql data) may not be fully pre - 
served in the Kepler pipeline processing (van Cleve 2010 ). 
That is not to say that all long-term trends are removed 
from the data by the Kepler processing, but the vari- 
ability statistics provided by NStED (and used in this 
study) are more sensitive to variability shorter than a few 
weeks. The primary effect of the Kepler pipeline is over- 
correction for shorter datasets (like the QO data) and 
fainter stars, but the pipeline is also capable o f adding or 
enhancing variability within the light curves (van Cleve 
20T0| 



Because we are interested in the overall variability 
statistics of the sample and not in the variability or pe- 
riodicity of any one individual star, the sheer size of the 
sample (~ 150, 000 stars) helps alleviate the specific ef- 
fects of any one star. In addition, the variability statis- 
tics presented in this work are in reasonable ag ree ment 
with s t atisti cs presented by Jenkins et al. (20101 and 



van 



Cleve (2010), and also in reasonable agreeme nt with the 
variability statistics of jBasri et al. ( 2010a|b l, who use a 
"range" of variability to describe the statistics. However, 
the results presented here should be viewed as a prelim- 
inary exploration of the public data set and are subject 
to revision as the Kepler project matures and improves 
the data products. 

2.2. Sample Segregation 

To help understand the variability statistics, we have 
utilized the Kepler Input C atalog (KIC 

2010) to separate 



et al. 2005 Batalha et al. 



stars into broad spectral and 



Latham 
tEe 



ummosity classes. The 
KIC includes stellar parameters (temperature and sur- 
face gravity) derived from photometric observations 
{u,g,r,i,z,DD051,J,H,Ks); a "Kepler Magnitude" 
corresponding to the bandpass of the instr ument is de- 
rived from the ground-based photometry (Koch et al. 
20101. The primary purpose of the KIC was to identity 
FTTj, and K (and M) dwarfs and separate them from the 
background giants in the field by utilizing photometry 
to determine line-of-sight extinction, effective tempera- 

''' http://nsted.ipac.caltech.edu 
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tures, and surface gravities (see |Batalha et al." (2010) for 
a description of the KIC algorithms and target selection 
process) . These derived values are available as part of the 
KIC information attached to each Kepler time series file. 
Of the 152,919 light curves available, 143,221 stars have 
KIC temperatures and surface gravities which we have 
used to separate the sample into dwarfs and giants by 
surface gravity and into spectral classes by temperature. 

The KIC temperatures and surface gravities are based 
upon isochrone fittin g utilizing the ATLAS9 models 



([Batalha etU.] |2010[ ). The KIC survey utiUzed the 
DD051 hiter which is sensitive to the MgH+Mgb line 
strength whic h varies as a function o f surface gravity for 



G and K stars ( [Majewski et al.|2000| . |Basri et al.| ( |2010b I 

showed that the KIC did a reasonably good job of sepa- 
rating giants from dwarfs, particularly for the G and K 
stars which dominate the sample. 

Separating the dwarfs and giants with a single value of 
surface gravity was not found to be sufficient. For exam- 
ple, a single surface gravity cut at \og{g) = 4.0 produces a 
bimodal distribution of the surface gravities for the giant 
star distribution and a truncated tail for the dwarf dis- 
tribution of surface gravities; these artificial structures in 
the distributions indicated that the giant sample was sig- 
nificantly contaminated by dwarf stars at the 20% level. 
In an effort to transition more naturally between giants 
and dwarfs, we have employed a three-section (empirical) 
surface gravity cut determined from the surface gravity- 
effective temperature HR diagram (see Figure IT]). For 
three separate temperature ranges, a star was consid- 
ered to be a dwarf if the surface gravity was greater than 
the value specified in the following algorithm: 



log(g) > 



(2.8 X io-'*r, 



eff) 



if Teff > 6000 
if Teff < 4250 
if 4250 < T^ff < 6000 



The delineation between dwarfs and giants is shown in 
Fig. IT] by the dashed line with the dwarfs and giants 
highlighted in blue and red, respectively. 

The total number of stars separated into dwarfs and gi- 
ants are 126,092 and 17,129, respectively. There is clear 
separation in the distributions of surface gravity for the 
two groups of stars (see middle panel Fig.[l|. The median 
surface gravities for the dwarfs and giants are, respec- 
tively, log((?) — 4.5 and log(g) = 3.0 with a small overlap 
in surface gravity near log{g) — 3.7. The overlap is likely 
dominated by sub-giants but represents a small contam- 
ination rate for both the dwarf and giant samples. The 
temperature distributions of the stars are shown in the 
lower panel of Fig. [l] the median dwarf and giant temper- 
atures are 5500 K and 4800 K, respectively. The dwarfs 
and giants have further been separated into temperature 
bins correspon ding roughly to the spectral types A, F , 



We have specifically explored the contamination rate of 
the M dwarfs with giant stars, by placing the M dwarfs on 
a 2MASS JHKs color-color diagram, where the dwarf and 
giant colors are sufficiently different to enable separation 
(see Figure [4|. Note that all of the M-dwarfs, as identi- 
fied from the KIC, have surface gravities of \og{g) > 4; 
yet, it is clear from the color-color diagram that a frac- 
tion of those identified as dwarfs are indeed giants. Using 
J — H = 0.75 mag as the boundary between dwarfs and 
giants, we find that only « 4% (108/2460) of the en- 
tire sample of stars identified as M-dwarfs in the KIC 
actually have infrared colors of a giant star. However, 
these contaminating stars are overwhelmingly brighter 
than the general M-dwarf sample with 80% of the giant- 
color "dwarfs" having a Kepmag brighter than 13.5 mag 
(see Fig. |4]). Thus, at the bright-end of the M-dwarf 
sample (Kepmag< 13.5 mag), the giant contamination 
rate is > 50% (87/170). The inverse contamination is 
also evident. The entire M-giant sample is much smaller 
with only 23 stars in total, but, of these, 6 (~25%) have 
JHK colors of dwarfs. The contaminating M dwarfs are 
systematically fainter than the true M-giants. For the 
sake of uniformity and continuity, we have not moved 
the contaminating sources into corresponding "correct" 
category; we do, however, exclude them when calculating 
the variability fractions. 

3. VARIABILITY 

For the spectral and luminosity classes defined above, 
we have assessed the distributions of the dispersion and 
variability to understand the broad stellar variability 
characteristics across the stellar spectrum. The analysis 
presented here utilized the statistics provided by NStED 
where the data were assessed using the native 30-minute 
sampling and the full 33-day time baseline of the dataset. 
The time series data are characterized by the dispersion 
about the median (am) and by the reduced chi-square 
assuming a constant median value for the light curve 
(x^). The first part of the study discusses the measured 
dispersions, and the second part of the study discusses 
the variability fractions of the stars within each group of 
stars. We also briefiy explore the variability fraction of 
the stars as a function of the time baseline of the dataset 
(1-33 days). 

3.1. Photometric Dispersion 

Figure [5] shows the 30-minute, 33-day photometric dis- 
persion as a function of Kepler magnitude for all the stars 
and separated out by dwarfs and giants, and Figure|6]dis- 
plays the dispersions to the same scale, but separated by 
temperature as well. The grey dashed lines in Figs.[5]and 
[6] correspond to the upper boundary on the uncertainties 
determined empiricall y for a constant back ground com- 
ponent (see Cupper 



Jenkins et al. 20101. The grey 



G, K, and M (|Johnson||1966| |Drilling & Landolt||2000|); solid line represents the median uncertainty value as a 

function of Kepler magnitude determined from the un- 
certainties provided with the data product as part of the 
light curves. To give some quantitative context to the 
numbers of stars within Figs. [5] and [6) Table [2] tabulates 
the numbers of stars within four different ranges of dis- 
persion. As instrumental precision plays a key role in the 
dispersion of the stars (particularly at the faint end) , the 
stars are grouped not only by stellar class but also by 
magnitude range. 



the temperature binning for each spectral class is listed 
in Table [T] and illustrated in Figure [2j The numbers 
of A and M stars are relatively small in comparison to 
the F, G, and K stars, but are maintained in the study 
for completeness. The temperature distributions clearly 
show that the G and K stars (and F dwarfs) dominate 
the sample. The magnitude distributions of the stars, 
separated by temperature and by dwarfs and giants, are 
shown in Figure |3] 
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There are a few specific aspects to tlie dispersion di- 
agrams that are worth noting. At fainter magnitudes 
(Kepmag > 14 mag), the model uncertainties (grey sohd 
and dashed hnes in Figures [51 [6]) track the stellar dis- 
persio n distribution fairly well (see also Jenkins et al. 
20101. At brighter magnitudes (Kepmag 14 mag), the 



model uncertainties track the lower bound of the mea- 
sured dispersions suggesting that Kepler is nearing the 
noise floor of the stars. This effect is most clearly seen 
in the giant stars where the Kepler data have sufhcient 
precision to detect the floor of the variability for the gi- 
ant stars. The G- and K-giants occupy a very narrow 
range of photometric dispersion between 0.1 — 1.0 mmag 
- completely independent of the magnitude. This nar- 
row range of dispersion is most clearly apparent in the 
dispersion distribution histograms (Figure^. 
The ubiquity of variability in g iants has been noted 



previously (Gilliland et al. 



2008) for a set of galactic 



bulge stars obse r ved b y HS'l over a time span of 7 days 



Gilliland et al. (20081 found the typical amplitudes of 
variability was ~ 0.5 mmag for the G giants and in- 
creased to ^ 3.5 mmag for the late-K to early-M giants. 
We see a very similar trend in the dispersion which is 
most clearly demonstrated in Figure [8] where we have 
plotted the photometric dispersion as a function of ef- 
fective temperature. While there is a scattering of stars 
with large dispersions (and the number of M giants is 
very small), the giant stars occupy a very narrow region 
of variability that is correlated with temperature. As 
expected from stellar ev olution, the larger, cooler g iants 
are more variable (e.g., Kjeldsen & Bedding 19951 and 
the variability spans two orders of magnitude (0.1 — 10 
mmag) . 

The dwarf stars are more complicated to interpret be- 
cause their intrinsic dispersion is on the order of (or less 
than?) the photometric precision. Taken as a whole, they 
are more quiescent than t he giant stars, as expect ed and 



demonstrat e d previously (Gilliland et al. 2008; Jenkins 
et al.||2010l Ivan Cleve||20l0| iBasri et al.||2010b|). But 



|Basri et ai:||2010b^ 
,t all magnitudes (Fie 



there is a sample of stars at all magnitudes (l*'ig. [5| 
and all temperatures (Figs. [6j [7j |8| where the average 
dispersion is ^ 5 mmag. Histograms of the dispersion 
(Fig. [t]) and plotting the dispersion as a function of tem- 
perature (Fig. [8| highlight the bi-modal dispersion, but 
show that only a relatively small percentage of stars are 
in the higher dispersion region. 

Visual inspection of a sample of 50 light curves (10 
light curves per temperature bin) in the high dispersion 
region indicates that these light curves are often peri- 
odic. Utilizing the NStED online periodogram servic^ 
~ 90% of the inspected light curves displayed one or 
more significant periods (the origin and distribution of 
the periods were not explored in this work). A similar 
visual inspection of 50 stars in the lower dispersion re- 
gion (but flagged as variable with xl > 2) revealed that 
the variability was dominated by more stochastic "white 
noise" rather than periodic variability, and only ~ 25% 
of the stars displayed signiflcant periodicity. It is also 
possible that these higher dispersions are an artifact of 
the data processing; however, this bi-modal dispersion 
distribution (Fig. | 7| for the dwarfs is also reported (al- 
though weaker) in Basri et al. (2010b) where they report 



a "variability excess" for those stars t hat are periodic ver- 
sus those stars that are not periodic (Basri et al. (2010b I 
independently processed the data utilizing an empirical 
polynomial fltting process). The details of the variability 
(e.g., periodic or stochastic, amplitude, structure) have 
not been fully explored in the work presented here, but 
it should be noted that variability do es not necessarily 
mean periodic behavior ( Howell|[2008 ) as all the stars in 
the visual inspection were flagged as variable, but not all 
stars were (obviously) periodic. 

The Kepler light curves are precise enough that even 
small variations in the light curve can lead to high dis- 
persion, where in typical ground-based transit survey 
data, the dispersion would remain relatively unchanged; 
transiting Jupiter-sized planetary companions can signif- 
icantly affect the measured dispersion. To help put this 
into perspective, we have over-plotted the positions of 
the known Kepler-fleld plane ts (BOKS-1, Hat-P7, TrES- 
2, Kepler-4,5,6,7,8 ; Howell e raL][20T0l |Pal et al.| [2008l 



' http:/ /nsted. ipac.caltech.edu/applications/ETSS/kepler_index 



O'Donovan et al.||2006, Borucki et al.||2010[ ) on the dis- 
persion diagrams (i'lgs. 5j8 ) . The dispersions of the light 
curves in these systems are ^ 2 mmag, except for Kepler- 
4 where the light curve dispersion is ~ 0.2 mmag. These 
light curves are nearly flat, to within the noise, except 
for the deep exoplanetary transits. If the transits are re- 
moved and the light curve statistics are recalculated, the 
dispersions decrease by almost an order magnitude for all 
the light curves except Kepler-4. All the planets (except 
Kepler-4) are Jupiter-sized with transit depths of ~ 1%), 
and it is the transits which dominate the statistics of the 
light curves. Kepler-4 is a much smaller (Neptune-sized) 
planet with a transit depth of only ^ 0.1% which is com- 
parable to the overall dispersion of the light curve. 

3.2. Variability Fractions 

The photometric dispersions (33-day baseline, 30- 
minute sampling) alone are not sufficient to assess the 
fraction of stars that are variable as the dispersion is de- 
pendent on the apparent magnitude of the targets, and, 
in particular, the dispersion for the dwarfs is at (or near) 
the precision limits of the instrument (for the 30-minute 
cadence). A more natural statistic is the reduced chi- 
square (xt) which takes into account the uncertainties 
(as reported in the public data light curve files). This 
analysis makes use of the provided point-to-point un- 
certainties in the light curves. For all light curves, the 
reduced chi-squares are calculated for the 33-day base- 
line (30-minute sampling) with respect to a constant me- 
dian value and are plotted as a function of temperature 
(Fig. |9]). For variability assessment purposes, a star is 
considered just-barely variable if > 2, significantly 
variable if xt > 10: ^^id very variable if xl > 100. A 
~ 2 corresponds to an excess dispersion of approxi- 
mately 1.5 times that of the measurement uncertainties; 
a « 10 corresponds to an excess dispersion of approx- 
imately 3 times that of the measurement uncertainties, 
and a x^ ~ 100 corresponds to an excess dispersion of 
approximately 10 times that of the measurement uncer- 
tainties. 

The measured fractions of stars that are variable are 
dependent upon the brightnesses of the stars as the 
instrumental precision decreases as the stars become 
fainter. Figure 10 plots the variability fractions as a 
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function of the Kepler magnitude for each of the stel- 
lar subgroups. The magnitude bins are 0.5 magnitudes 
in width, and the fractions were calculated for the three 
reduced chi-square categories listed above. Uncertain- 
ties on the fraction s were calculate d using standard error 
propagation (Everett et al. 2002 1 . Some of the magni- 
tude bins have very few stars particularly at the bright 
end, and this is reflected in the relatively large error bars. 
At the very bright end of the Kepler sample (Kepmag 

< 11 mag), the variability fractions for the stars with 

> 2 are all near unity indicating that the Kepler 
precision at 30-minute sampling is approaching the 33- 
day noise floor for the stars. Also for the brightest stars, 
the fractions of stars that are significantly variable {x1 > 
10) is 50 — 100% depending on the sub-group of stars. 
The fractions decrease as the stellar magnitude increases; 
this is, of course, a direct result of the decrease in the 
instrumental precision as the stars become fainter. Not 
surprisingly, at all brightnesses, the fractions of stars that 
are significantly variable {xt > 10—100) are less than the 
fraction of stars that are just-barely variable {x1 > 2). 
For the dwarfs, as the stars grow fainter (Kepmag > 14 
mag), the variability fractions are typically dominated 
by the extremely variable stars {xt > 10 — 100), again a 
result of the lower precision on the fainter stars. 

A summary of the variability fractions is given in Ta- 
ble [3l where the fractions for each group of stars have 
been calculated for the entire sample (Kepmag < 16 
mag) and for the brighter end of the sample (Kepmag 

< 14 mag). The table includes variability fractions for 
the whole light curve baseline (33-days) as well as for 
1-day and 10-day time baselines which are discussed in 
more detail in lXO 

3.2.1. The Dwarf Stars 

The G-dwarfs are the least variable group of dwarf 
stars with > 80% of the stars being stable (Kepmag < 16 
mag; xt < 2); even with the magnitude restricted to the 
brightest stars (Kepmag < 14 mag), the variability frac- 
tion of the G-dwarfs is ~ 30%. The floor for the G-dwarf 
dispersion appears near 0.04 mmag (see Figl6|. The K- 
dwarf and F-dwarf stars have comparable variability frac- 
tions with « 50% of the stars identified as a variable if 
the magnitude is restricted to Kepmag < 14 mag. The 
F-dwarfs have a higher variability fraction if F-dwarfs of 
all magnitudes are considered, but this is a result of the 
different magnitude distributions of the K and F dwarfs 
in the sample (see Fig. |3| . The relative number of F and 
K stars brighter than and fainter than Kepmag ~ 14 
mag differ, with the K-dwarfs having significantly more 
fainter stars than brighter stars. The differing magni- 
tude distributions for the F- and K-dwarfs is a result of 
the target selection criteria optimized for searching for 
transiting planets around as m any appropriate stars as 
possible. dBatalha et"aI1|2010l ) . 

The M-dwarts, not surprismgly, are less stable than 
the G, K, and F dwarfs. If the magnitude is restricted to 
Kepmag < 14 mag, nearly 70% of the stars are variable 
(this is after the removal of stars which appear to have 
giant star colors of J — H > 0.75 mag). The variability 
fraction drops to 36% if all the stars in the sample are 
considered. The A-dwarfs have a similar variability frac- 
tion with sa 70% of the A-dwarfs being variable. Similar 



results are seen in the Hipparcos variability statistics, 
where the A-dwarfs display the highest vari ability frac- 
tions (see Figure 2 of [Eyer fc Mowlavi||2008"| . The large 
fraction of variable A-stars is likely the result of the A- 
star group (as identified from the KIC) including stars in 
the instability regime, such 7 Dor, 5 Scu, slowly pulsat- 
ing B-stars (SPBs), RR Lyr, and /3 Cep stars. Hippar- 
cos found the variability fractions of these sub-groups to 
range from 10 - 100%. 

3.2.2. The Giant Stars 

At the precision of Kepler, nearly all of the giants are 
variable, with 94%, 99%, and 100% variability fractions 
(33-day baseline, 30-minute sampling) for the G, K, M 
giants, respectively. The 6 M-stars with dwarf J-H colors 
have been removed from the statistics. The G-giants 
variability fraction is slightly reduced by the faint end 
of the brightness distribution, where the stability floor 
approaches the instrument limit for Kepmag « 14 mag, 
but only because the stability floor of the G giants is 
0.1 mmag versus 0.3 mmag for the K-giants. For the M- 
giants, the dispersion and stability floor is substantially 
higher at levels of ~ 10 mmag. 

The variability fraction of the giants f ound in the Ke- 
pler data i s consistent with the work of Gilliland et al.| 
pOOSj) and |Eyer fc Mowlavi] ([2008 ) , where the majority of 
the giants were found to be variable and a strong correla- 
tion of variability with decreasing temperature along the 
gia.nt bra nch was found. In ground-based work (Henry et 
al.|2000 ), a similar trend was found, but the photometry 
was not precise enough (^ 1 mmag) to see the variability 
of the hotter G- and early-K giants (< 0.5 mmag). The 
timescales of the variations in these works were found to 
be inconsistent with rotational modulation of a spotted 
photosphere, and were found to be more consistent with 
acoustic oscillations of the atmospheres, with the varia- 
tions of the late-K and M giants consistent with radial 
pulsations, and the variations of the more stable G and 
early-K giants dominated by non-radial pulsations. 

Assuming that the photometric dispersion in the Ke- 
pler giants is also dominated by acoustic oscillations, the 
photometric variations can be used to predict radial ve- 
locity amplitudes of the oscillations. Kjeldsen & Bedding 
( |1995 !) developed a calibrated relationship between the 
velocity of oscillations and the photometric amplitude 
variations: 



20.1 X 10- 



0.55Aimy V5777; 



where arv is the oscillation velocity of the star, {/S.F/F)\ 
is the photometric flux change at the observed wave- 
length A, and Tg// is the effective temperature of the 
star. Using this relation, we have calculated the expected 
radial velocity oscillations for the G- and K-giants based 
upon their photometric dispersions and effective temper- 
atures (see Figure 11). 

The bulk of predicted radial velocity dispersions are 
centered around 10 — 20 m/s with 90% of the velocities 
< 30 m/s. The K-giants have a symmetric distribution 
centered at {(Jrv) « 20 ± 5 m/s. This is in go od agree- 
m ent wi th a radial velocity study of K-giants (Frink et 
al. 12001 ) where it was found that the radial distribution 
of K-giants could be described with a Gaussian of mean 
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20 m/s and width of 11 m/s with a long tail to higher 
velocity dispersions. The agreement with the predicted 
and measured distributions for representative samples of 
K-giants suggests that the variability observed by Kepler 
is dominated by acoustic oscillations in the atmospheres 
of the giants. 

The G-giants predicted velocities show a bi-modal 
structure with peaks near 10 and 20 m/s with the 
stronger peak towards lower radial velocity variations. 
The magnitude distributions of the G-giants that have 
predicted radial velocity amplitudes of < 15 m/s and 
those that have predicated radial velocity amplitudes of 
> 15 m/s are indistinguishable indicating that the bi- 
modality is not related to the brightness (and hence, the 
photometric precision) of the stars, but rather is intrinsic 
to the sample. The radial velocity appears uncorrelated 
with temperature, but does appear to have a weak anti- 
correlation with surface gravit5Q suggesting that the G- 
giant sample may contain a sampling of dwarfs and sub- 
giants, which are atmospherically more stable than the 
G-giants. 

3.2.3. Time Dependent Variability 

The analysis, thus far, has been performed on the 
full 33.5 day time baseline of the quarter-1 dataset (30- 
minute cadence), but in reality, stars are variable on 
a variety of timescales depending on the source of the 
variability (e.g., flares , pulsations, rotation, and eclipses; 
Eyer fc Mowla.vi||2008[ ). A full detailed study of the vari- 
ability as a function of the time baseline and sampling 
rate is beyond the intent and scope of this paper, but we 
have briefly explored how the variability fractions change 
depending upon the length of the dataset investigated. 
It should be noted (as discussed above) that the Kepler 
public product may remove long- term variability or en- 
hance some forms of variability 
the detailed results of this short st 
with that in mind. 

For each of the light curves, we have assessed the light 
curve properties with progressively longer time baselines 
starting at 1 day and extending to 33 days. The me- 
dian, the dispersion about the median, and the reduced 
chi-square assuming a constant median value were calcu- 
lated for each time interval. To help alleviate biases that 
might arise from sampling the light curves in progres- 
sively longer samples in a single direction, the statistics 
were calculated by sampling the light curves in the time- 
forward direction (0—1 day, — 2 day — 3 day, . . . 
— 33 day) and in the time-backward direction (33 — 32 
day, 33 — 31 day 33 — 30 day, ... 33-0 day) and the 
results were averag ed. In Table \3\ and Fi gure |12| the 
dependency of the derived variability fractions are sum- 
marized. As with the overall variability statistics, the 
analysis was performed for all stars (Kepmag < 16 mag) 
and for the brighter end of the sample (Kepmag < 14 
mag). 

The overall fractions of giants that are variable {xt > 
2) do not change to within the uncertainties of the frac- 
tions as the time baseline is increased from 1 day to 33 

^ The Kendall-r non-parametric rank correlation value between 
the surface gravities and the predicted radial velocity oscillations is 
—0.75 (number of standard deviations from zero is 100); a value 
of —1 would indicate a perfect anti-correlation. 
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days. The fraction of giant stars that are more signifl- 
cantly variable {xl > 10 and xl > 100) does grow by 
~ 4 — 5% from the 1-day baseline to the 33-day baseline. 
The small growth of the variability fractions is likely a 
result of the fact that nearly all of the giants are ob- 
served to be variable at the precision of Kepler and the 
variability fraction has little room to change. 

For the dwarf stars, the overall variability fractions 
{xl > 2) increase by « 1—5%, as the baseline is increased 
to 33 days. As with the giants, the variability fraction 
changes more substantially for those stars that are more 
signiflcantly variable {xl > 10 and xt > 100). Larger 
amplitude variability requiring longer time periods is not 
surprising and ha s been observed previously (e.g., Eyer 



fc Mowlavl] 2008 ) . Increasing the time baseline from 1 
day to 33 days increases the variability fractions for the 
M-dwarf, K-dwarf, G-dwarf stars variability more than 
for the F-dwarf and A-dwarf stars. This suggests that 
the lower mass stars are predominately characterized by 
variability with timescales of weeks (e.g., rotational mod- 
ulation) while the higher mass stars are predominately 
characterized by variability with timescales of days (e.g., 
pulsations) . 

3.3. Galactic Distribution 

The Kepler field spans approximately 12 degrees in 
galactic latitude (6 « 8° — 20°). Over this range of lati- 
tude, the different galactic populations may play a role in 
the variability fractions. Because the target samples are 
mostly magnitude-limited, the differing intrinsic bright- 
nesses of the stars lead to differing median distances of 
the stars for each sub-group, and hence, to differing me- 
dian hei ghts (z) above the galac tic plane for a given line 
of sight. Walkowicz et al. (2010) found a higher fraction 
of the flaring M and K dwarfs at lower z- heights and they 
suggested that they were sampling primarily the young 
thin disk. Their work inspired us to try to understand 
the overall variability fraction of the sample as a function 
of latitude and z- height for each of the stellar sub-groups. 

A subset of the Kepler Field was selected (see Figure 
13) to remove the effects of the rotation of the Kepler 
Held with respect to the Galactic plane. The median 
temperature and magnitude for each category of stars 
was used to determine a "typical" distance for the stars, 
assuming zero attenuation by interstellar dust (see Table 
The z- height of each star was computed from the typ- 
ical distance for its sub-group, its apparent magnitude, 
and its galactic latitude. This simple estimation assumes 
that each star within a subgroup has the same absolute 
magnitude. While this, of course, is not strictly correct, 
the typical spread of absolute magnitude within a sub- 
group is 1 — 2 mag, corresponding to only a factor of 
1.2—1.5 in the distance. The z-height distributions of the 
stars (Figure 14 ) mostly follow the expected exponential 
decay for a disk of the form N oc exp(—z/zo) where Zp 
is the characteristic scale height of the disk (Ciardi et 
al.||1996l lJuric et al. 2008). Each z-height distribution 



was fitted with a decaying exponential and the resulting 
scale heights are listed in Table |4j The exponential fits 
work best for the intrinsically brightest stars (e.g., the 
F-dwarfs, A-dwarfs, K-giants and G-giants) where local 
distribution effects are minimized. 

The M-dwarfs and K-dwarfs, with distances of only a 
few hundred parsecs and scale heights of z < 100 pc. 
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are dominated by stars located nearer to the disk plane 
and by stars within the solar neighborhood. The G-, 
F-, and A-dwarfs all display larger characteristic scale 
heights (z = 100 — 180 pc), but are all within the ex- 
pected size of the young thin disk. The K- and G-giants 
have scale heights oi z = 200 — 250pc which is charac- 
teristic of the older thin disk (Juric et al. 2008). If the 
stars came from only this one disk population, it is ex- 
pected that ^ 90% of the stars will have z-heights within 
z < 2.3zo. The actual fractions are listed in Table [4j all 
of which are significantly below 90%, indicating that the 
thick disk may contribute to the overall sample - partic- 
ularly at higher galactic latitudes. The thick disk has a 
scale height of ~ 900 pc and a scaling fraction of ^ 10% 
( [Juric et al.l|2008[ ). 

it the thin disk contributes only a portion (albeit the 
majority fraction) to the sample of stars observed by Ke- 
pler, a variation in the variability fraction as a function 
of galactic latitude (i.e., scale height) might be expected. 
Figure [l5] displays the fraction of stable stars (x^ < 2) 
and variable stars (x^ > 2) as a function of galactic lat- 
itude for each of the sub-groups (K-giants and M-giants 
are not included in this sample as "all" of the K-giants 
are variable at the precision of Kepler). The M, K, G, 
and F dwarfs all show an increase in the variability frac- 
tion as the galactic latitude gets lower (i.e., closer to 
the plane). Moving higher in galactic latitude, the vari- 
ability fractions decrease by ^ 0.1% over the 10° span 
of the Kepler Field. This could indeed be the result of 
sampling younger stars in the plane at lower latitudes as 



young stars are expected to be more active (West et al. 
20081. Indeed, the flaring rate of M-dwarfs as a lunc- 



tion of z-height suggests that stars located nearer to the 
galactic plane are more active and, hence, more variable 
(Walkowicz et al. 20101, in reasonable agreement with 
what is discussed here. 

An alternative explanation is that the background con- 
tamination is higher when looking closer to along the 
galactic plane and that the increased variability is the 
result of more significant blending of the primary star 
with fainter background stars. The slope of the vari- 
ability fraction as a function of latitude is strongest for 
the low luminosity stars (M- and K-dwarfs), weakens as 
the intrinsic luminosity of the stars increases (G- and 
F-dwarfs), and is not apparent for the most intrinsi- 
cally bright stars (A-dwarfs and G-giants). As the Ke- 
pler sample is magnitude limited with similar magnitude 
ranges for each of the stellar sub-groups, the different 
sub-groups are essentially sampling different distances 
(see Tabled. 

For the thin disk (zo ^ 300 pc), the path length to 
outside the disk (z ~ 600 pc) is sa 4300 pc at 6 ^ 8°, but 
only « 1750 pc at 6 ~ 20°. For the M-dwarfs with typical 
distances of 200 pc, the latitude-change corresponds to a 
background path length (for a conic volume) difference of 
nearly 40% from low (& ~ 8°) to high (6 ~ 20°) galactic 
latitude. For the G-, F-, and A-dwarfs {d ^ 1000 pc), the 
background volume difference is < 20%. The reduction 
in background path length is approximately 50% from M- 
dwarfs to F-dwarfs, which is also the fraction by which 
the slopes of the variability fraction vs latitude change 
from M-dwarfs to G- and F-dwarfs (see Figure 15). The 
A-dwarfs do not display a reduction in the variability 



fraction at higher latitudes; if anything, they exhibit a 
weak (and somewhat insignificant) increase in variability 
at higher latitudes. The G-giant stars, with typical dis- 
tances that are larger than the line of sight distances to 
the "top" of the exponential disk at 6 ^ 10 — 20°, show 
no dependence of the variability fraction on the galactic 
latitude. All of this is consistent with background stars 
contributing to the variability of the primary stars. 

Without a full model of the stellar galactic distribution 
coupled with a priori knowledge of the true variability 
fraction of the relative populations, it is difficult to disen- 
tangle these scenarios (true variability fractional changes 
as a function of latitude vs. changes in the background 
contamination rate). However, t he apparent correlation 
of flare rates with lower z- height ( Walkowicz et al.|2010 l 
does suggest that the higher variability fraction at lower 
galactic latitudes may be real and the result of sampling 
a systematic younger population. 

3.4. M-Dwarf Variability 

M-dwarfs are favorable targets to search for earth-sized 
planets because the transits are relatively deep (~ 1 — 3 
mmag), and the radial velocity signatures are relatively 
large (~ 10 m/s). In addition, planets in the habitable 
zones of M-stars are in relatively short orbits (10 — 20 
days) compared to that of the habitable zones for sun-like 
stars (~ 1 year). As a result there has been a strong in- 
terest in the community for searching for planets around 
M-sta rs (Irwin et al. 2009; Charbonneau et al. 2009';'Beaii| 
et al.pOlO) . Thus, understanding the M-dwarf variabiF 
ity aniplitudes and fractions is critical to understanding 
how complete such transit and radial velocity surveys can 
be. 

In the previous sections (j 3.1|3.2 ), the overall variabil- 
ity fraction of the M-dwarfs was toimd to be ^ 40 — 70% 
with dispersions of f7„i ~ 3 — 5 mmag, depending on 
the brightness of the stars being considered. As an al- 
ternative, in this section we identify a small sample of 
relatively bright, certain M-dwarfs based on well- vetted 
proper motion catalogs and analyze their variability in 
more detail. These M dwarfs include all of the M dwarfs 
in the Kepler field (with Ql light curves ) from the Gliese 
and LHS catalogs ( Stauffer et al.|[20T0 ), and the bright- 
est stars in the LSl^M catalog with V — J > 2.6 (i.e. 
colors consistent with an M dwarf). A plot of J — vs 
H — K confirms that these are indeed M dwarfs (Fig- 
ure 16 ). Only four of these stars have KIC Tg/ / or \og{g) 
- the rest would be absent from statistical studies which 
rely on T^ff and \og{g) to identify dwarfs, and, of these, 
one would have been classified as a giant. 

To understand the variability of these bright M-dwarfs 
on the time scales relevant to planetary transits, we have 
calculated the short-term 12-hour variability for each of 
the light curves, by computing the dispersions in running 
12-hour time bins. The median of all the 12- hour bin 
dispersions for each light curve were calculated and taken 
as representative of the 12-hour variability timescales for 
the M-dwarfs. The dispersions for the full time series (33 
day) and for the 12-hour timescales are listed in Table [5j 
In all cases, the dispersion on the 12-hour timescale is 
smaller than the full 30 day dispersion, and for many of 
the stars, the dispersion drop s to the photometric limit 
of instrument (see Figure 17). 

On the 33-day timescaleptne dispersion is bimodal with 
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peaks near 0.1 mmag and 5 mmag. The 0.1 mmag peak is 
dominated by stars which are quiet to the precision of the 
instrument, and « 1/2 (15/29) of the sample are variable 
with > 2 and dispersions of am ^ 1 mmag. For the 
12-hour timescale, the variability fraction drops signifi- 
cantly with only 6 stars that have dispersions am > 0.5 
mmag. The high dispersion is likely caused by rotational 
variability with periods of 1 day or longer; thus, it is not 
all that surprising that the dispersion drops when the 
light curves are sampled at 12-hour timescales. A prime 
example of this is LHS6343 (KIC 100022 61) which is a 
newly discovered transiting brown dwarf 
20101. The dispersion for the entire ligh 



Johnson et al. 



curve IS 



mmag, but the 12-hour timescale dispersion matches the 
out-of-eclipse dispersion of « 0.7 mmag - much like what 
is obs erve d for the transiting planets around FGK stars 
(see 



3.1) 



4. SUMMARY 

An analysis of the variability statistics of the stars in 
the Quarter- 1 publicly released Kepler data has been per- 
formed. The Kepler data cover 33.5 days and are sam- 
pled at a 30 minute cadence. The Kepler Input Cata- 
log parameters have been used to separate the 150,000 
stars into dwarfs and giants which were further separated 
into temperature bins corresponding roughly to spectral 
classes A, F, G, K and M. 

The majority of the dwarf stars were found to be photo- 
metrically quiet down to the per-observation (30 minute) 
precision of the Kepler spacecraft. The derived variabil- 
ity fractions range from 10 - 100% depending on the stel- 
lar group and brightness range explored. The G-dwarfs 
are the most stable with < 20% of the all the stars in the 
sample having a > 2. The G-dwarfs appear to have a 
dispersion noise floor of ~ 0.04 mmag for the 30-minute 
sampling of the Kepler data. 

At the precision of Kepler, > 95% of K, G, and M gi- 
ants are variable with noise floors of ^ 0.1 mmag, ~ 0.3 
mmag, and ~ 10 mmag, respectively. The photometric 
dispersion of the giants is consistent with acoustic varia- 
tions of the photosphere. The photometrically-predicted 
radial velocity distribution for the K-giants is in agree- 
ment with the measured distribution; the G-giant ra- 
dial velocity distribution is bimodal which may indicate 
a transition from sub-giant to giant. 

We also briefly explored the dependence of the vari- 
ability fractions as a function of time baseline of the light 
curves. In general, increasing the length of the light curve 
baseline increased the fraction of stars that are variable. 
For the dwarf stars, the lower mass stars were found 
to be predominately characterized by variability with 
timescales of weeks (e.g., rotational modulation) while 
the higher mass stars were found to be predominately 
characterized by variability with timescales of days (e.g., 
pulsations). For the giant stars, the variability fractions 
changed very little from a 1-day sampling to a 33-day 
sampling. 

A study of the distribution of the variability as a func- 
tion of galactic latitude suggests sources closer to the 
galactic plane are more variable. The scale height dis- 
tribution of the dwarfs is consistent with the young thin 
disk, and the scale height of the giants is consistent with 
the older thin disk. For the lower mass stars (M, K, and 
G dwarfs) , the variability fraction decreases with increas- 



ing galactic latitude. This may be the result of sampling 
differing populations as a function of latitude and prefer- 
entially sampling younger stars at lower galactic latitudes 
within the Kepler field. 

In addition to the statistical study of M dwarf vari- 
ability using the 2500 relatively anonymous probable 
M dwarfs in the Kepler field, we have also examined 
the variability of 29 known M dwarfs in the Kepler field 
drawn from the GJ, LHS, and LSPM catalogs. The anal- 
ysis of the known M dwarfs indicates that the M dwarfs 
are primarily variable on timescales of weeks presum- 
ably dominated by spots, rotation, and binarity. But 
on shorter timescales of hours-to-days, the stars are qui- 
eter by nearly an order of magnitude. At these shorter 
timescales, the variability fraction of the M-dwarfs drops 
from ^ 40% to ^ 20%. The shorter timescales are rel- 
evant for searches of planetary transits which typically 
last a few hours. In general, a search for transiting earth- 
sized planets around M-stars should not be hampered by 
the typical stellar variability of M-dwarfs. 

The authors would like to acknowledge the referee 
for his or her extremely insightful and useful comments 
which made this a better paper. Portions of this work 
were performed at the California Institute of Technol- 
ogy under contract with the National Aeronautics and 
Space Administration. This research has made use of 
the NASA/IPAC Star and Exoplanet Database, which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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Table 1 

KIC-Based Temperature Bins 



Spectral 


Dwarf 


Dwarf 


Giant 


Giant 


Type 


Te//Range 


Number'' 


Tj,yy Range 


Number'* 


A 


> 7300 


2311 (2296) 







F 


(6000 - 7300 


23750 (15996) 







G 


(5300 - 6000 


66682 (17940) 


> 4800 


9880 (9877) 


K 


(4000 - 5300 


30889 (4874) 


(3800,4800] 


7226 (7225) 


M 


< 4000 


2460 (171) 


< 3800 


23 (17) 



'^ Number is parentheses is number of stars brighter than Kepmag< 14 mag. 



Table 2 

Stars Separated by Class, Magnitude, and Dispersion" 



Stellar 


Kepler 


# of Stars 


# of Stars 


# of Stars 


# of Stars 


# of Stars 


Group 


Magnitude 


in Magnitude 


with 


with 


with 


with 




Range 


Range 


0- < 0.1 


0- = 0.1 - 1 


0- = 1 - 10 


0- > 10 


M Dwarfs'" 


<10 


3 








3 







10-12 


13 


2 


3 


8 







12-14 


154 





83 


57 


14 




14-16 


2182 





1503 


529 


150 


K Dwarfs 


<10 


18 


3 


10 


3 


2 




10-12 


264 


75 


83 


94 


11 




12-14 


4588 


92 


3063 


1212 


221 




14-16 


26019 


1 


19892 


5184 


942 


G Dwarfs 


<10 


63 


26 


27 


10 







10-12 


1297 


716 


257 


294 


27 




12-14 


16566 


542 


13376 


2332 


316 




14-16 


48756 





43533 


4550 


673 


F Dwarfs 


<10 


141 


28 


81 


29 


3 




10-12 


1948 


490 


966 


429 


54 




12-14 


13906 


402 


11407 


1806 


291 




14-16 


7755 





7158 


519 


78 


A Dwarfs 


<10 


231 


114 


56 


56 


5 




10-12 


739 


280 


210 


226 


20 




12-14 


1326 


124 


658 


460 


84 




14-16 


15 





5 


8 


2 


M Giants'^ 


<10 


6 








3 


3 




10-12 


8 








4 


4 




12-14 


3 








1 


2 




14-16 

















K Giants 


<10 


350 





273 


73 


1 




10-12 


1683 


1 


1468 


200 


11 




12-14 


5192 


1 


4776 


349 


66 




14-16 

















G Giants 


<10 


233 


5 


209 


17 


2 




10-12 


1619 


38 


1467 


93 


21 




12-14 


8025 


7 


7689 


255 


74 




14-16 


















^ Dispersions in milli-magnitudcs [mmag]. 

^ The 108 contaminating giants in the M-dwarf sample have been removed from the statistics. 
The 6 contaminating dwarfs in the M-giant sample have been removed from the statistics. 



Table 3 

Variability Fractions^ 



Time 








> 2 










> 10 






xl 


> 100 




scale 


Category 


< 16 


mag 


< 14 


mag 




< 16 


mag 


< \A 


mag 


< 16 


mag 


< 14 


mag 


33 day 


All Dwarfs 


0.269 


(0.002) 


0.461 


(0.004) 





177 


(0.001) 


0.269 


(0.003) 


0.123 


(0.001) 


0.200 


(0.002) 




M Dwarfs'' 


0.367 


(0.014) 


0.690 


(0.083) 





285 


(0.012) 


0.497 


(0.065) 


0.209 


(0.010) 


0.474 


(0.064) 




K Dwarfs 


0.298 


(0.004) 


0.532 


(0.013) 





224 


(0.003) 


0.347 


(0.010) 


0.161 


(0.002) 


0.316 


(0.009) 




G Dwarfs 


0.183 


(0.002) 


0.325 


(0.005) 





125 


(0.002) 


0.199 


(0.004) 


0.090 


(0.001) 


0.164 


(0.003) 




F Dwarfs 


0.421 


(0.005) 


0.555 


(0.007) 





212 


(0.003) 


0.279 


(0.005) 


0.127 


(0.003) 


0.169 


(0.004) 




A Dwarfs 


0.705 


(0.023) 


0.705 


(0.023) 





559 


(0.019) 


0.558 


(0.019) 


0.434 


(0.016) 


0.435 


(0.016) 




All Giants 


0.962 


(0.011) 


0.962 


(0.011) 





717 


(0.009) 


0.717 


(0.008) 


0.210 


(0.004) 


0.210 


(0.004) 




M Giants': 


1.000 


(0.343) 


1.000 


(0.343) 


1 


000 


(0.343) 


1.000 


(0.343) 


1.000 


(0.343) 


1.000 


(0.343) 




K Giants 


0.996 


(0.017) 


0.996 


(0.017) 





886 


(0.015) 


0.886 


(0.015) 


0.314 


(0.008) 


0.314 


(0.008) 




G Giants 


0.938 


(0.014) 


0.938 


(0.014) 





593 


0.010) 


0.593 


(0.010) 


0.135 


(0.004) 


0.135 


0.004) 


10 day 


All Dwarfs 


0.268 


(0.002) 


0.461 


(0.004) 





167 


(0.001) 


0.260 


(0.003) 


0.097 


(0.001) 


0.167 


(0.002) 




M Dwarfs'' 


0.367 


(0.014) 


0.678 


(0.082) 





269 


(0.012) 


0.474 


(0.063) 


0.163 


(0.009) 


0.404 


(0.058) 




K Dwarfs 


0.299 


(0.004) 


0.535 


(0.013) 





212 


(0.003) 


0.341 


(0.010) 


0.127 


(0.002) 


0.275 


(0.008) 




G Dwarfs 


0.182 


(0.002) 


0.328 


(0.005) 





117 


(0.001) 


0.189 


(0.004) 


0.069 


(0.001) 


0.128 


(0.003) 




F Dwarfs 


0.417 


(0.005) 


0.551 


(0.007) 





205 


(0.003) 


0.272 


(0.005) 


0.102 


(0.002) 


0.141 


(0.003) 




A Dwarfs 


0.690 


(0.022) 


0.689 


(0.023) 





548 


(0.019) 


0.548 


(0.019) 


0.417 


(0.016) 


0.417 


(0.016) 




All Giants 


0.961 


(0.011) 


0.962 


(0.011) 





713 


(0.009) 


0.713 


(0.008) 


0.202 


(0.004) 


0.202 


(0.004) 




M Giants': 


1.000 


(0.343) 


1.000 


(0.343) 


1 


000 


(0.343) 


1.000 


(0.343) 


1.000 


(0.343) 


1.000 


(0.343) 




K Giants 


0.995 


(0.017) 


0.995 


(0.017) 





883 


0.015) 


0.838 


(0.015) 


0.303 


(0.007) 


0.304 


'0.007) 




G Giants 


0.937 


0.014) 


0.937 


0.014) 





589 


0.010) 


0.589 


0.010) 


0.128 


0.004) 


0.128 


0.004) 


1 day 


All Dwarfs 


0.264 


(0.002) 


0.432 


(0.004) 





089 


(0.001) 


0.183 


(0.002) 


0.042 


(0.001) 


0.097 


(0.002) 




M Dwarfs'' 


0.320 


(0.013) 


0.567 


(0.072) 





128 


(0.008) 


0.281 


(0.046) 


0.077 


(0.006) 


0.152 


(0.032) 




K Dwarfs 


0.279 


(0.003) 


0.491 


(0.012) 





082 


(0.002) 


0.205 


(0.007) 


0.033 


(0.001) 


0.093 


(0.005) 




G Dwarfs 


0.193 


(0.002) 


0.312 


(0.005) 





054 


(0.001) 


0.105 


(0.003) 


0.021 


(0.001) 


0.050 


(0.002) 




F Dwarfs 


0.397 


(0.005) 


0.513 


(0.007) 





152 


(0.003) 


0.211 


(0.004) 


0.077 


(0.002) 


0.108 


(0.003) 




A Dwarfs 


0.678 


(0.022) 


0.679 


(0.022) 





519 


(0.018) 


0.520 


(0.019) 


0.389 


(0.015) 


0.340 


(0.015) 




All Giants 


0.955 


(0.010) 


0.955 


(0.010) 





668 


(0.008) 


0.668 


(0.008) 


0.174 


(0.003) 


0.174 


(0.003) 




M Giants': 


1.000 


(0.343) 


1.000 


(0.343) 


1 


000 


(0.343) 


1.000 


(0.343) 


0.941 


(0.323) 


0.941 


(0.328) 




K Giants 


0.995 


(0.017) 


0.995 


(0.017) 





839 


(0.015) 


0.838 


(0.015) 


0.266 


(0.007) 


0.266 


(0.007) 




G Giants 


0.925 


(0.013) 


0.925 


(0.013) 





544 


(0.009) 


0.544 


(0.009) 


0.107 


(0.003) 


0.107 


(0.003) 



Values in parentheses are uncertainties based upon the propagation of errors of the counting statistics. 
'' The 108 contaminating giants in the M-dwarf sample have been removed from the statistics. 
': The 6 contaminating dwarfs in the M-giant sample have been removed from the statistics. 



Table 4 



Galactic Distributions 




Median 


Median 


Typical 


Scale 


Fraction of 




Tcff 


Kepmag 


Distance 


Height Zo 


Stars 


Category 


[K] 


[mag] 


[pc] 


[pc]'' 


z < 2.3zo 


M Dwarfs 


3800 


15.3 


200 


35 


0.97 


K Dwarfs 


5000 


15.1 


600 


75 


0.68 


G Dwarfs 


5700 


14.7 


1000 


105 


0.56 


F Dwarfs 


6200 


13.6 


1100 


185 


0.83 


A Dwarfs 


8000 


12.3 


1200 


165 


0.77 


G Giants 


5000 


13.1 


2800 


235 


0.40 


K Giants 


4700 


12.7 


2500 


215 


0.43 



^ Calculated by fitting the z-height distributions in Figure |14[ 
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Table 5 

M-Dwarf Stars 



Star 
Name 


KIC 
ID 


KIC 

[K] 


KIC 

log(,g) 

[cm s 


Cm 

33 day 
[iiimag] 


12 hour 
[mmag] 


LHbbool 


2104791 






6.09 


1.69 


LbF1912 + OOZD 


o o o n o A 






0.42 


0.41 


T O T) 1 rir\r» i ori i n 

LoFlyUy+oylU 


4U4jcioy 


3713 


4.385 


7.94 


0.17 


OrJ4U99 


A 1 A on 1 o 

414291o 






4.05 


0.10 


G J4113 


44 ^U9o / 






0.06 


0.04 


LbF1917+4UU7 


OUUiociD 






0.10 


0.10 


LbF 1947+4020 








2.93 


0.13 


TCDinoc 1 At 1 f\ 

LbF19oo+4119 


6049470 






2.05 


0.09 


T CTDiniA 1 /IIT? 

LoF 1919-1-4127 


D117dU2 






3.84 


2.13 


T oni o O 1 A1 AT 

LoFlo5o+4147 


Do45oo5 






2.84 


0.08 


LbF195o+4149 


ATI O O CT 


3201 


0.07 


0.20 


0.18 


LbF19z7+4zol 


7Uoc)D7U 




■ ■ ■ 


0.36 


0.28 


LbF1944+4zo2 


7U494DO 


4033 


4.505 


1.47 


0.08 


LbF191z+4zo9 


Ti n/i; on'7 

71UboU7 






0.12 


0.11 


ijoFiyiz+4olO 


/oyo9iu 






0.20 


0.19 


T uoc o /in 

LllbDo49 


7siU5J5 






0.34 


0.33 


T0T>10Cr/1 1 A A A'7 

LbF 1854+4447 


SDU7728 






7.06 


0.31 


LSP2001+4500 


8846163 






0.17 


0.15 


LHS3429 


8872565 






0.12 


0.11 


LSP1933+4515 


8957023 


3553 


4.117 


7.82 


0.18 


LHS3420 


9201463 






39.8 


7.13 


GJ1243 


9726699 






11.6 


9.69 


LHS6343 


10002261 






3.07 


0.68 


LSP1857+4720 


10258179 






0.11 


0.10 


LSP1854+4736 


10453314 






0.17 


0.14 


GJ4083 


10647081 






4.69 


0.08 


LSP1916+4949 


11707868 






0.11 


0.14 


LSP1948+5015 


11925804 






0.79 


0.55 


LSP1919+5130 


12555642 






1.92 


0.10 
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Figure 1. Top: KIC-based Surface Gravity - Effective Temperature HR-diagram of the stars in the analysis sample. The dashed black 
line marks the delineation to separate dwarfs (blue) and giants (red). Center: Histograms of the surface gravity for the dwarfs (blue) and 
giants (red). The vertical dashed lines mark the median surface gravity values. Bottom: Histograms of the effective temperatures for the 
dwarfs (blue) and giants (red). The vertical dashed lines mark the median temperature values. 
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Figure 2. Temperature distributions (binsizc = 100 K) of the stars selected as dwarfs (fop) and giants {bottom). The color-coding 
illustrates the separation of the dwarfs and giants into temperature groups (i.e., spectral types). 
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6 8 10 12 14 16 18 6 8 10 12 14 16 18 

Kepmag [m^g] Kepmag [mag] 





Kepmag [mag] Kepmag [mag] 

Figure 3. Kepler magnitude distributions of the stars in the sample. Dwarfs and giants are represented by the blue and red hashed 
histograms, respectively. The panels represent the different temperature groups as labeled in the figures. 
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Figure 4. Top: 2MASS color-color diagram for the stars identified as M-dwarfs (left) and M-giants (right) based upon the KIC surface 
gravities and effective temperatures. The green-hashed area marks the main sequence; the blue-hashed area marks the giant branch, and the 
red-hashed area marks the L-dwarf locus. The diagonal lines mark the reddening zone for typical galactic interstellar extinction (R=3.1). 
Bottom: Magnitude distributions for stars identified as dwarfs (left) and as giants (right) by their surface gravity. The black histograms are 
for stars with dwarf-like J-H colors; the red histograms are for stars with giant-like J-H colors. These plots show that the M-stars brighter 
~Kepmag< 13.5 mag are predominately giants, regardless of their KIC classification. 
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Figure 5. Photometric dispersion (30 minute sampling; 33 day timescale) of each star is plotted as a function of magnitude for all stars 
(top), for just the dwarfs (center), and for just the giants (bottom). In the top and center plots, the locations of the 7 known planets in 
the sample are shown (red: BOKS-1, Hat-P7, and TrES-2; green: Kepler-4,5,6,7,8). The grey line re presents the rnedian uncertainty as 
reported in the Kepler data product. The dashed grey curve is the uncertainty upper limit curve from jjenkins et al.| ( |20l"oj l. 
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Figure 6. Photometric dispersion (30 minute sampling; 33 day timescale) of each star is plotted as a function of magnitude separated out 
by temperature and surface gravity as labeled in each panel The solid grey line represe nts the median uncer tainty value as reported in the 
Kepler data product. The dashed grey curve is the uncertainty upper limit curve from |Jenkins et al.| l |20l"o| . 
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Figure 7. Distributions of tlie (logaritlimic) photometric dispersion (30 minute sampling; 33 day timescale) separated out by effective 
temperature and luminosity class as labeled in each panel. 
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Figure 8. Photometric dispersion (30 minute sampling; 33 day timescale) of each star is plotted as a function of effective temperature 
separated out by temperature (colors and labels) and surface gravity (top and bottom panels). The black points in the top panel mark the 
locations of the seven known planets in the sample (+: BOKS-1, HAT-P7, and TrES-2; X: Kepler-4,5,6,7,8). 
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Figure 9. Photometric reduced chi-square (30 minute sampling; 33 day timescale) of each star is plotted as a function of effective 
temperature separated out by temperature (colors and labels) and surface gravity (top and bottom panels). The black points in the top 
panel mark the locations of the seven known planets in the sample (+: BOKS-1, HAT-P7, and TrES-2; X: Kepler-4,5,6,7,8). 
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Figure 10. Variability fractions of stars as a function of tiie brigiitness (Kepler magnitude). The contaminating giants in the M-dwarf 
sample and the contaminating dwarfs in the M-giant sample have been removed from the statistics. The blue curves represent the fractions 
of stars with > 2; the green curves represent the fractions of stars with xi > 10; the red curves represent the fractions of stars with 
Xl > 100. 
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Figure 11. Distribu tion of radial velocity osci llations of G (blue) and K (red) giants predicted from the photometric dispersion and 
effective te mperat ure jKjeldsen &: Bedding|1995| l. The dashed line marks the median radial velocity oscillation for the K-giant sample of 
[Prink et al.]pOOl| . 
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Figure 12. Variability fraction distributions as a function of sampling. Each panel represents a different group of stars. There are two 
panels for each group; one panel for all the stars in the sample, and one panel where the stars were restricted to a Kepler magnitude of 
14 or brighter. The contaminating giants in the M-dwarf sample and the contaminating dwarfs in the M-giant sample have been removed 
from the statistics. The blue curves represent the fractions of stars with xi > 2; the green curves represent the fractions of stars with 
xi > 10; the red curves represent the fractions of stars with x2 > 100. 
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Figure 13. Galactic coordinates plot of the positions of all the stars in the sample. The red-box delineates the 10° X 10° region used to 
explore the variability as a function of galactic latitude. 
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z [pc] z [pc] 

Figure 14. z-height distributions for tlie dwarfs and G and K giants. The black smooth curves represent the best fit exponential curves 
to the distributions. The M-giants have been excluded from the plot because of the low number (23) in the sample. 
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Figure 15. Galactic latitude distributions (binsize=l°) for the dwarfs and G giants. The black curves represent the fraction of stars 
within that galactic latitude bin that are deemed "stable" (xj < 2), and the red curves represent the fraction of those stars that are deemed 
"variable" {xi > 2). The black dashed line is a best fit to the variability fraction as a function of galactic latitude with the parameters of 
the line fit given in each panel. 
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Figure 17. Left panels show the photometric dispersion plotted as a function of magnitude for the KIC-identified M-dwarfs with a 
color restriction of J-H< 0.75 mag (black dots) and for the outside identified M-dwarfs (red). The solid grey line represents the median 
uncertainty as reported in the Kepler data product. The dashed grey curve is the uncertainty upper limit curve from Jenkins e t al., | ,2010| . 
Right panels show the distributions of the (logarithmic) photometric dispersion (binsize = 0.2 dex) for the known M-dwarts (red points in 
left figures). The top panels reflect the dispersion of the known M-dwarfs determined for the entire light curve (30 days); the bottom panels 
reflect the dispersion calculated from the point-to-point differences on 12-hour timescales (only for the known (red points) M-dwarfs). The 
KIC-identified M-dwarfs (black dots) are shown at the 30-minute cadence dispersion in both plots for reference. 



